Vortex quantum tunnelling versus thermal activation in ultrathin superconducting 

nanoislands 



W. V. Pogosovi^2 and V. R. Misko^ 
^ Departement Fysica, Universiteit Antwerpen, Groenenborgerlaan 171, B-2020 Antwerpen, Belgium and 
^Institute for Theoretical and Applied Electrodynamics, 
Russian Academy of Sciences, Izhorskaya 13, 1254-12, Moscow, Russia 

We consider two possible mechanisms for single- vortex fluctuative entry /exit through the surface 
barrier in ultrathin superconducting disk-shaped nanoislands made of Pb and consisting of just 
few monoatomic layers, which can be fabricated using modern techniques. We estimate tunnelling 
probabilities and establish criteria for the crossover between these two mechanisms depending on 
magnetic field and system sizes. For the case of vortex entry, quantum tunnelling dominates on the 
major part of the temperature/flux phase diagram. For the case of vortex exit, thermal activation 
turns out to be more probable. This nontrivial result is due to the subtle balance between the 
barrier height and width, which determine rates of the thermal activation and quantum tunnelling, 
respectively. 
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Recent progress in miniaturization techniques has 
made it possible to create ultrathin quasi-two- 
dimensional superconducting structures, which consist of 
just few monoatomic layers, and even one-layer films. In 
contrast to previously studied mesoscopic superconduc- 
tors both thermal and quantum fluctuations can be 
expected to be important in such extremely thin struc- 
tures [^l. In addition, strong confinement of condensate 
in lateral dimensions leads to a rather rich variety of pos- 
sible behaviors. In Ref. Q, vortices in superconducting 
nanoislands of Pb with thickness of few nanometers have 
been studied experimentally. It was observed that hys- 
teresis effects are significantly weaker than expected from 
the theory. In another experiment 0], similar islands 
have been explored. The experiment has revealed a full 
suppression of the surface barrier for the sirwle-voitex 
entry/exit. These results (see also Refs. d, Q) suggest 
that some kind of fluctuations of the superconducting 
state might be responsible for the observed features. In- 
deed, it was show in Ref. [81 that the surface barrier 
height, under the conditions of the experiment Q, is low 
enough to be suppressed by thermal fluctuations. 

The aim of the present paper is to consider another 
possible mechanism of surface barrier suppression in ul- 
trathin superconducting nanoislands, which is due to 
quantum fluctuations, as well as to establish a domi- 
nant mechanism, i.e., thermal activation (TA) or quan- 
tum tunnelling (QT). We consider a disk-shaped super- 
conducting island consisting of just few monoatomic lay- 
ers with radius much larger than the zero-temperature 
coherence length. The magnetic fleld is applied perpen- 
dicular to the island surface. Presence of the boundaries 
results in the surface barrier which prevents vortex entry 
or exit. In the absence of fluctuations, the barrier for 
vortex entry disappears, when the magnetic fleld reaches 
some critical value. The barrier for the vortex exit is sup- 
pressed at a lower magnetic fleld. The difference between 
these two flelds thus reflects the signiflcance of hysteresis 
effects. In the London model, which remains applicable 



down to low temperatures, these two critical flelds can be 
calculated by using the approach developed by Buzdin 
P . Note that macroscopic QT of Abrikosov vortices has 
been studied both experimentally and theoretically for 
the case of long current-carrying strip in zero applied 
fleld, where vortex penetration is due to the current (lOl - 
IT^ . as well as for bulk superconductors, where energy 
barriers appear due to the pinning of vortices by spatial 
inhomogeneities [l3| or due to intrinsic pinning induced 
by the layered structure [ij] (for the QT of Josephson 
vortex, see recent papers [15|). In the present paper, we 
apply a combination of known approaches to the particu- 
lar case of the ultrathin mesoscopic disks, which became 
experimentally testable quite recently. Apart from the 
interest from the viewpoint of a possible experimental 
examination, it is also attracting to see how strong con- 
flnement in lateral dimensions affects an interplay be- 
tween quantum and thermal scenarios, since it is known 
that probabilities of these processes are controlled by two 
different parameters, which are barrier width and barrier 
height [13, [H, [il] , these parameters being strongly de- 
pendent on system sizes and geometry. 

In order to estimate the probability of macroscopic QT 
of Abrikosov vortex to/from the island, we use an ap- 
proach developed by Caldeira and Leggett 
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which 

allows one to take into account energy dissipation. For 
a vortex with mass rriy in a potential V{r) (to be found 
within the London model), the Euclidean action is given 

by 117] 
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The last term in the right-hand side (r.h.s.) of Eq. ([IJ, 
which is nonlocal, takes into account the energy dissipa- 
tion. In the reciprocal space, we obtain 
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It is convenient to introduce a frequency-dependent ef- 
fective mass as meff{uj) = rriy + 77/|a;|. Note that it is 
shown below that mass is much smaUer than the con- 
tribution from the viscosity. 

In principle, one can use Eq. ([2]) to find optimal u{uj) 
exactly. However, it is known that quite accurate results 
for the effective action can be obtained by using dimen- 
sional estimates 13] (in addition, we also use the method 
of Ref. ji3). Following Refs. [iS [H, we find the char- 
acteristic tunnelling frequency ujq from the balance of 
kinetic and potential energies as 
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where VJ, and lb are barrier height and width, respectively. 
Eq. is a quadratic equation, from which luq can be 
found. In the viscous limit {myVb ^ v'^lb)^ ^^e solution 
reads as 
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where we have kept a linear term. It will be seen below 
that the dissipative regime is indeed relevant to the sit- 
uation we consider here. Using Eq. ([2]) and (jH), we find 
the action as 
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Note that in the case of quadric-plus-cubic potential, 
which was solved exactly in Ref. |16| . there appears a 
prefactor 7r/2 instead of 2 in the dominant term of the 
r.h.s. of Eq. ([5]), which demonstrates a high accuracy of 
the method of dimensional estimates. 

The viscosity coefficient for the quasi-two-dimensional 
superconductor, in the case of Ohmic dissipation, is 
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where J? at is a two-dimensional resistance, which depends 
on island thickness d as[ll| Rn — Rodo/d, where dn is 
the thickness of one monolayer. According to Ref. 
Rq is nearly equal to 6 kOhm for Pb. 

The mass of a vortex can be estimated using the well- 
known expression by Suhl [l8l |. It includes two contri- 
butions: (i) due to the kinetic energy of the vortex core 
and (ii) due to the electromagnetic energy. In Ref. [l9| , 
another contribution was found, which was attributed 
to the shear deformations of the crystal lattice. This 
additional contribution is typically of the same order as 
those in Ref. 18| , so that it does not violate viscous- limit 
conditions. These conclusions apply to the dirty-limit 
regime, which is realized for superconducting nanostruc- 
tures placed on a disordered substrate, as relevant for the 
experiments In contrast, in the clean limit, one has 

to take into account a much larger vortex mass arising 
from the quantization of the electron states in the vor- 
tex core j20| , which may depend on the frequency of the 



external drive. In the dirty limit, however, it becomes of 
the order of Suhl mass [II [HI, [2|. 

The probability of QT, within exponential accuracy, 
is Pqt ^ exp{—SE/fi), while the probability of TA is 
PTA-expi-Vb/kBT). 

The effective action, in the dissipative limit, is mainly 
due to the width of the barrier, as follows from Eq. ^ 
(see also Refs. [13, [ll)- Thus, the probability of QT 
is controlled by barrier width. At the same time, the 
probability of TA is determined predominantly by the 
barrier height. Below we show that the subtle balance 
between the two characteristics of the barrier leads to a 
rather rich behavior of the system we here study. 

In the London model, the energy as a function of the 
vortex position p can be found analytically [9] . For a disk 
of radius R much smaller than the effective penetration 
depth Aj^ = X{T)^ /d, the position-dependent part of the 
energy reads as |9|: 
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where $ is a magnetic flux through the disk. 

We consider a superconducting island placed on a dis- 
ordered substrate 01 ■ Thus, the diffusive regime is re- 
alized with the mean free path nearly equal to 2d. It was 
shown experimentally in Ref. (23j that the Ginzburg- 
Landau theory remains applicable for such extremely 
thin films. For A(T) and f (T), we use the diffusive-limit 
expressions given by 
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where bulk Aq = 40 nm and ^0 = 80 nm for Pb. For the 
dependence of Tf, on d, we utilize the expression deduced 
from the experiments on few- monolayer Pb films (23l |: 
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where Tco = 7.2 K is the critical temperature of Pb in a 
bulk, while the " critical thickness" d^ corresponds nearly 
to 1.5 monolayers. 

It is easy to show from Eq. ^ that the position of the 
saddle point, which separates vortex-free and one- vortex 
states, is 
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The condition for the vortex entry can be obtained by 
equating to i? — ^(T). After simple algebra, we obtain 
the critical fiux ^in, at which barrier for vortex entry 
disappears 
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When magnetic flux is lower than the barrier height 
Uin for the vortex entry is nonzero. It can be calculated 
as the difference in energies in a saddle-point configura- 
tion and at p — R~ ^(T): 
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The condition for vortex exit is pc — 0, which im- 
mediately leads to the critical flux $out = 'i'o- When 
^ > ^out, the barrier height Uout for vortex exit is 
nonzero. It can be found as a difference of energies in 
a saddle-point state and at p = 0: 
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Thus, it is possible to find analytically barrier heights 
both for vortex entry and exit. However, it is not possible 
to find the barrier widths analytically, since it amounts 
solving transcendental equations. Namely, the barrier 
width for vortex entry is given by the difference in two 
values of p, the first one being R — £, (T) , while the second 
one corresponding to another solution of equation U{p) = 
U{R — £,(T)). Similarly, the barrier width for vortex exit 
is given by a nonzero value of p, at which U{p) — U{0). 
The above two equations are then solved numerically. 

Let us focus on the case of an island with i? = 120 nm 
and d = 1 nm. These are realistic parameters which 
are available in modern experiments. The disk radius in 
terms of ^(0) is 10, while d is nearly 3 monolayers. This 
leads to Tc = 3.6 K. In Fig. 1, we show typical energy 
profiles for vortex entry /exit. It is quite remarkable that 
the width of the barrier for the entry is generally much 
smaller than that for the exit. Since the probability of 
QT is controlled by the square of the barrier width, we 
may expect that such an asymmetry can provoke thermal 
escape of a vortex from the island to be more probable 
than QT. 

Now we are going to address two issues. First, we 
would like to understand if TA and/or QT is realistic for 
the system we consider. Second, we want to establish, 
which mechanism is dominant depending on T and $. 

Our calculations show that exponents for both QT and 
TA turn out to be ^ 10 — 100 almost for the whole phase 
diagram. We believe that values close to ten or even sev- 
eral times larger are sufficient for the vortex to overcome 
the barrier: Indeed, it was shown in Refs. [1, [13, HI] by 
evaluating the pre-exponential factors that similar values 
of the exponent result in observable rates for TA for both 
bulk and nanosized superconductors. In addition, we en- 
sure that we are indeed in the dissipative limit, since 
rriyVb/rf'll is several orders of magnitude smaller than 1. 

Our results are summarized in Fig. 2. Curve 1 shows 
$m as a function of T . Curve 2 represents $ojit- Dashed 
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FIG. 1: (Color online) Typical curves for the energy of the 
system as a function of the vortex position in a supercon- 
ducting island of radius 120 nm and thickness of 1 nm. Curve 
1 describes the barrier for vortex entry at $ = 2.1$o and 
T — O.GTc. Curve 2 corresponds to the barrier for vortex exit 
at $ = 1.5^0 and T — O.GTc- The energies are normalized by 
their values at saddle points. Dotted lines are guides for eyes 
used to depict barrier widths. 
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FIG. 2: (Color online) The "phase diagram" "Magnetic field 
vs. Temperature" for the island (see the text). The inset 
shows the dependence of the decimal logarithm of Se/^ on 
T/Tc along curve 4. 



curve 3 gives $, at which the energies of the state with 
vortex at the disk center and at p = R — f (T) are equal. 
For the domain on the phase diagram, which is between 
curves 3 and 1, vortex entry is energetically favorable. 
For the region between curves 3 and 2, vortex exit is 
preferable. 

Curve 4 separates two domains: below this curve 



4 



change of the topological charge is more probable through 
thermal fluctuations, while above this curve quantum 
fluctuations dominate. We denote the corresponding 
temperature as Tq, while /c_bTo ~ (^/2??)(Vf,/^fc)- For 
the domain between curves 4 and 1, QT of the vortex 
into the disk is more probable than TA. We see that the 
width of this region expands with decreasing T, which is 
reasonable. For the region confined between curves 1, 3, 
and 4, TA of vortex into the disk is more probable. It 
is remarkable that this scenario survives down to rather 
low T. The domain between curves 2 and 3 corresponds 
to the vortex escape from the island. In contrast to the 
case of vortex exit, it is fully determined by the ther- 
mal mechanism, even at relatively low T . This is due 
to the fact that barrier width for exit is quite large and 
is controlled mostly by the disk radius, while the one 
for the entry is mainly determined by ^(T) (see Fig. 1). 
The inset in Fig. 2 shows the dependence of Vh/ZcsT 
on TjTc along curve 4, where it is equal to SeI^- It is 
also of interest that curve 4 is actually oriented nearly 
horizontally in this "phase diagram". It crosses curve 1, 
which corresponds to the threshold for the vortex entry 
in the absence of fluctuations, at T w O.SSTc. This im- 
plies that, along curve 1, QT remains more probable up 
to quite high T. 

Finally, two dotted lines in Fig. 2 depict borders of the 
regions, close to the threshold values of where the ex- 
ponent for the dominant mechanism of vortex entry / exit 
is lower than 10, so that a change in the system's topo- 
logical charge becomes highly probable. These lines can 
be therefore interpreted as the ones corresponding to the 
disappearance of the barrier. They intersect in the point 
on the curve 3, where both barriers are suppressed, so 
that the vortex exit and entrance must become nearly 
reversible. 

In the present paper, we used the method of dimen- 
sional estimates, based on manipulations with two main 
parameters of the barrier, namely, barrier width and 
height. In order to provide an additional support for 
our results, we also applied the approach of Ref. [3], 
which was initially developed to address quantum creep 
in layered high- Tc materials. Within this method, the ef- 
fective action is first calculated in the high-temperature 
regime, which is purely due to TA. Then, the crossover 
temperature Tq is found by solving the equation of mo- 
tion taking into account only the first Fourier harmonic 
as a correction to the static solution. For T < Tq, the ac- 
tion is estimated, again, by using the TA expression but 
at r = Tq. We found that the results of the method used 
here turn out to be quite similar to the results of this ap- 
proach. For instance, Tq can be obtained from Eq. (25) 
of Ref. by dropping the irrelevant elastic term, as 
ksTo ='-hU"{pc)/2mj. Since U"{pc) ~ -Vb/ll, we ar- 
rive at the similar result for Tq. By evaluating U"{pc), 
one can also easily show that fc^To ~ ($ — $o)) in agree- 
ment with curve 4 in Fig. 2. The difference in numerical 
prefactors of Tq, as found by two approaches, is essen- 
tially compensated by a very sharp slope of Tq as a func- 



tion of 

The important question is, how the predicted features 
can be tested experimentally. It is possible to detect 
presence of a vortex in a disk by using STM methods, 
as already has been done in Ref. y, |j] for smaller and 
thicker disks. In particular, Ref. Q studied the differ- 
ence between magnetic fields for which vortex entry and 
exit occur. It was shown that this quantity is signifi- 
cantly smaller than it could be expected from the theory 
when disregarding fluctuations. In our point of view, it 
is attractive to address the temperature dependence for 
average fields of vortex entry and exit. One may assume 
that the barrier becomes suppressed when the exponent 
approaches some critical value, the reasonable estimate 
being 10. As seen from Fig. 2, in the absence of fluc- 
tuations, the flux of vortex exit must be independent 
on T, while in the presence of fluctuations it must grow 
with increasing T. At low T, the corresponding value of 
magnetic flux must be close to $ojit- By expanding the 
r.h.s. of Eqs. (HO]), ([12]), and (HH) in powers of $ - ^out, 
one can easily show that the width of the fluctuative 
region where thermal fluctuations are dominant, varies 
with temperature as ^ VT. The same applies to the 
case of vortex entry. Contrarily, when quantum fluctu- 
ations are dominant, the width of the fluctuative region 
is nearly T-independent (except for a rather weak de- 
pendence through ^(T) and A(T)), so that it tends to a 
constant at T = 0. Hence, by analyzing the low-T behav- 
ior, one might distinguish between two different mecha- 
nisms of barrier suppression through fluctuations. For 
the vortex exit, we predict a crossover between QT and 
TA (see the upper dotted curve in Fig. 2). In this case, 
the slopes of the dependencies of the effective critical flux 
on T, which are nearly linear, must be different for low 
and high T. 

Perhaps, a more elaborate method to study TA of a 
single vortex in ultrathin islands has been proposed very 
recently in Ref. f26| , where it was concluded that islands 
of this kind provide an ideal platform to study such phe- 
nomena. It was suggested to change the applied field pe- 
riodically and slowly and to measure not only the average 
fields for vortex entry/exit, but also their distributions. 
The distribution profile must broaden with the increase 
of T in the case of TA. This feature provides an additional 
tool to distinguish between TA and QT scenarios. 

In summary, we have studied the possibility of both 
the macroscopic quantum tunnelling and thermal activa- 
tion of Abrikosov vortex through the surface barrier in 
ultrathin superconducting islands made of Pb and con- 
sisting of just few monoatomic layers. Such islands can 
be fabricated and studied using modern techniques. We 
have found that the barrier for vortex entry, as well as 
for the exit, can be suppressed in this system, even not 
necessarily in the immediate vicinity of threshold val- 
ues of magnetic flux in the absence of fluctuations. The 
dominant mechanism of vortex fluctuative entry/exit is 
determined by the interplay between the barrier height 
and width, which are the main quantities defining prob- 
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abilities of thermal activation and quantum tunnelling in 
a dissipative limit, respectively. A "phase diagram" of 
the island in the plane of magnetic flux and temperature 
was constructed. 
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